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Abstract 

Inertial magnetization dynamics have been predicted at ultrahigh speeds, or frequencies ap¬ 
proaching the energy relaxation scale of electrons, in ferromagnetic metals. Here we identify inertial 
terms to magnetization dynamics in thin Ni 79 Fe 2 i and Co films near room temperature. Effective 
magnetic fields measured in high-frequency ferromagnetic resonance (115-345 GHz) show an addi¬ 
tional stiffening term which is quadratic in frequency and ~ 80 mT at the high frequency limit of 
our experiment. Our results extend understanding of magnetization dynamics at sub-picosecond 
time scales. 
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The magnetization M{t) in ferromagnetic materials is generally understood to evo 
without memory of its prior motion. As described by the Landau-Lifshitz (LL) equation 
magnetization dynamics dM{t)/dt can be written in terms of the magnetization M{t) alone, 
excluding any temporal derivatives M(t)/df^. Magnetization is then ’inertialess’: it re¬ 
sponds to a step magnetic held H{t) with an instantaneous change in speed and with inhnite 
acceleration. 

As pointed out by Ciornai et al. and Fahnle et a/.0-4|, followed by other theoreticians j^, 




the absence of inertia is questionable for magnetization dynamics at very high frequen¬ 
cies. The high frequency behavior, > 100 GHz, becomes relevant in ultrafast (< ps-scale) 


magnetization switching schemes envisaged for magnetic information storage technology 




ll[ |. In metallic ferromagnets, the electrons which comprise the magnetization themselves 


possess inertia, expressed through a hnite ^loch-state) lifetime r and cannot change their 
momenta inhnitely quickly. These authorshave proposed that the LL equation should 
be amended to include an inertial term: 


dxw y. 

—— = -/io7[m X Heff am x (m x Heff)] ± am x t—— 
dt dt^ 


( 1 ) 


Here m is the unit magnetization vector, 7 is the gyromagnetic ratio, Hgff is the ef¬ 
fective held, and a is the Gilbert damping coefhcient. The additional nonlinear term 
±am X Td^m./d?t is second-order in time. It is straightforward to show that this inertial 
term leads to an effective held = =FQ:i^^'7‘/7, quadratic in frequency, which acts on 

magnetization dynamics. Increasing the frequency limit for ferromagnetic resonance (FMR) 
experiments by a factor of 10 will therefore increase the magnitude of the inertial ehective 
held by a factor of 100 . 

The sign of the inertial term in Eq. (1) refers to two diherent origins for the proposed 


inertial dynamics. Giornai et al. 


12 | argue that additional energy terms for the precession. 


typically neglected, come about from transverse magnetization motion. Finite rotational 
moments of inertia about the two transverse axes were introduced semiclassically to account 


for this energy in Ref. [^. On the other hand, Fahnle et a/.[3|, 1^ have pointed out that 
similar dynamics arise from a more complete consideration of damping. Here the inertia 
can be associated with the linear momentum of electrons: electrons cannot change their 
momenta k more rapidly during precession than their Bloch state relaxation time tb, and 
tend to remain at higher energy states. In this sense, inertia is proportional to tb- The 
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rotational momeYA of inertia (”+” in Eq. 1) softens the resonance frequency while the linear 
inertia due to relaxation in Eq. 1) stiffens the resonance frequency. 

In this work, we test for the existence of inertial effective helds which act on the mag¬ 
netization of the device ferromagnets permalloy (Ni 79 Fe 2 i, or Py) and Co. To do so, we 
measure FMR up to 345 GHz, roughly an order of magnitude higher than is accessible 
using conventional sources. Through measurements of the frequency-dependent resonance 
field fioHresi^^), we hnd in all samples a negative nonlinear shift quadratic in frequency 
which reflects a positive effective held added by the inertial term. The sign of the shift 
shows the dominance of inertia due to the hnite Bloch state relaxation timejs, 4| over the 
inertia due to hnite rotational moment of inertia of electrons a. nail samples the ex^acted 
values of r are comparable with low-huence limits of remagnetization times te measured in 
ultrafast demagnetization experiments. Our results show that high-frequency FMR (HF- 
FMR) provides a near-ground-state (~ 1 meV) view of the processes involved in ultrafast 
magnetization dynamics, heretofore accessible through far-from-equilibrium techniques. 


Py and Co thin hlms were deposited on Si/Si02 substrates using ultra-high vacuum mag¬ 
netron sputtering. The capping layers are Cu(5 nm)/Si02(5 nm) for Py and Cu(5 nm)/Ta(3 
nm) for Co; the seed layers are Ta(5 nm)/Cu(5 nm) in both series. The thicknesses tpM = 
6 , 10 and 30 nm for Py and 6, 10, 15, 30 nm for Co, respectively. Two series of held-swept, 
variable-frequency FMR experiments were carried out, both with applied held and static 
magnetization perpendicular to him plane. At low helds and frequencies (3-26 GHz, 1-2.5 
T, 298 K), a coplanar waveguide was used to measure the power transmission through the 
thin hlms. These measurements were carried out at Columbia University using an Fe-core 
electromagnet 13]. At high helds and frequencies (115, 230 and 345 GHz, 4-13 T, 273-280 
K) continuous-wave high-held electron paramagnetic resonance was used to measure the 
power rehection from the samples. These measurements were carried out at the LNCMI-G 

n 

using a superconducting magnet |141|. Six held sweeps were taken for HF-FMR to reduce 
measurement and hysteresis errors. 

Several steps have been taken to ensure the consistency of low- and high-frequency mea¬ 
surements. These are essential to detect the nonlinear term in magnetization dynamics. 
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FIG. 1. (a) Low-frequency resonance fields for Py and Co < 26 GHz. Lines and dashed lines are 

fits to (u/y = fj.o{Hres — ^eff)- The upper panel shows the htting residuals of the lower panel, 
(b) High-frequency FMR absorption spectra for Co 15 nm. The (horizontal) held axis shows the 
difference of measured high-frequency resonance helds with expected values based on extrapolation 
of the low-frequency resonance helds without inertial effects. The solid data points indicate the 
resonance held fXoHres, with hts shown in the dashed curve. 


First, the magnetic fields o: 
using paramagnetic markers 


■ both low- and high-frequency FMR were calibrated in-situ 
15[ |. An accuracy of 0.01% is obtained, corresponding to ~ 1 


mT at 345 GHz (~ 13 T). Second, for low-frequency FMR, each sample normal has been 
independently aligned along the applied held by rotating the waveguide with respect to two 
axes (with < 0.2° precision) to maximize the value of fioHres at 3 GHz. The samples were 
aligned optically, to lower precision, in HF-FMR. Third, we have corrected for the small (~ 
20 K) temperature differences between the two apparati, through variable-temperature FMR 
measurements at Columbia University, which allow us to correct for temperature-dependent 
magnetization p,oTfe//(T). See the Supplemental Information for details. 

Fig. 1(a) shows the low-frequency fioHres as a function of a;/27r for Py and Co. Contin¬ 
uous and dashed lines are hts to the linear Kittel equation cu/y = /io(ihres — Tfe//)- The 
top panel shows the htting residuals, which are randomly distributed within ±1 mT. In 
order to examine whether resonance held ohset exists at high frequencies, we plot in Fig. 
1(b) the HF-FMR lineshapes (average of six measurements) for Co 15 nm, as a function of 
/io(hf — with the resonance held extrapolated from low-frequency FMR. 


4 











The solid circles indicate the resonance field ji^Hres extracted from the FMR lineshapes for 
the three frequencies. The dashed lines are a fit of the resonance positions to a quadratic 
nonlinear term to frequency. A negative curvature is found in the dashed lines. This implies 
a negative sign in Eq. (1) with positive effective inertial field. Similar effects are found in 
other Py and Co films. 
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FIG. 2. Deviation of /xoi^res from the linear Kittel equation for (a) Py and (b) Co. Main panels: 
HoHres\exp — (Aio-^e// + ^/'y)\fit solid data points; in curved lines. 7 and r 

are the fit parameters from Eq. (2). Upper panels: fitting residuals of main panels as a function 
of frequency. Insets: full plots of noHres and fits to Eq. (2). 


To quantify the nonlinear term from Eq. (1), we fit the resonance helds of both low- and 
high-frequency FMR measurements together into the linear Kittel equation, plus a quadratic 
term: 

HoHres = hoMe// u/j - aoUj\/^ (2) 

taking M^ff, 7 and r as ht parameters. The parameter r governs the magnitude of the 
inertial term. The Gilbert damping coefficients Oq are extracted by htting the low-frequency 
FMR linewidths into /ioAifi /2 = ^qAHq + 2aou/'j, where hqAHq is the inhomogeneous 
linewidth. To estimate errorbars for HF-FMR, we take the scatter of all six measurements 
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into the data sequence. 

Figure 2 contains the central result of this manuscript. In the main panels we show 
experimental resonance helds ^oHres and hts to Eq. (2) as a function of frequency u)/27r 
for (a) Py and (b) Co. The conventional effective held terms + u/'^ have been 

subtracted away; plotting the data as /UoFfreslexp — (ho^e// + ^h)\fit (solid data points) 
and aow^r/yl/it (curves) isolates the inertial contribution. The ht values are tabulated in 
Table 1; error bars denote Sa of uncertainty. We hnd a reasonable ht to the data through 
the introduction of the quadratic term. For both Py and Co, the values of r stay roughly 
at the same level between 10-30 nm: ~ 0.1 ps for Py and ~ 0.3 ps for Co. For Py 06 nm 
negligible shifts in /ioFfres is found, as will be discussed. 


Sample 

hoMe// 

9eff 

ao 

T(po-Rres) 

r(poARi/2) 

Py 06 nm 

0.915 T 

2.106 

0.0075 

*0.28±0.06 ps 

0.09±0.09 ps 

Py 10 nm 

0.935 T 

2.100 

0.0072 

0.13±0.03 ps 

0.11±0.04 ps 

Py 30 nm 

0.993 T 

2.101 

0.0069 

0.12±0.05 ps 

0.06±0.08 ps 

Co 06 nm 

1.409 T 

2.147 

0.0067 

0.47±0.08 ps 

0.19±0.07 ps 

Co 10 nm 

1.547 T 

2.153 

0.0053 

0.26±0.09 ps 

0.14±0.07 ps 

Co 15 nm 

1.596 T 

2.154 

0.0051 

0.32±0.03 ps 

0.15±0.05 ps 

Co 30 nm 

1.672 T 

2.155 

0.0056 

0.26±0.05 ps 

0.21±0.09 ps 


TABLE I. Fit parameters of experimental resonance fields ^QH^esi^) to the inertial term Eq. (2) 
and linewidths fioAHi /2 to Eq. (3) for Py and Co. Values of low-frequency Gilbert damping 
oo are found from linear fits of low-frequency linewidths. Errorbars of r represent 3a deviations. 
*Misaligment analysis has been applied; see text for details. 


A primary source of uncertainty may come from sample misalignment in HF-FMR measure¬ 
ments. Given the experimental geometry, the samples cannot be aligned independently of 
mounting. Canting the sample oh normal alignment in HF-FMR will also lead to reduced 
resonance helds as in Fig. 1(b). Misalignment will be manifested as a linear deviation of 
ftQHres\exp — {l^o^eff + ^h)\fit from being zero at low frequencies (a;/27r < 30 GHz) in 
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Fig. 2. For Co samples, the deviations at low frequencies are negligible, indicating that 
Co samples are well aligned perpendicular to biasing field. For Py the same is true, with 
the exception of the 6 nm sample. This sample shows a significant linear deviation at low 
frequencies (green circles), indicating that the sample plane may not be perfectly perpendic¬ 
ular to the biasing held. According to simulations, shown in the Supplemental Information, 
such a deviation would be produced by a misalignment of ~ 6° for Py 06 nm, which will 
add an additional negative fioHres shift of ~ 20 mT. With this correction, the estimated r 
is increased from 0.01 ps to 0.28 ps in Py 06 nm, consistent with the larger value for Co 6 
nm. The corrections are much smaller for other samples and have been neglected. 

The inhuence of the inertial dynamics which we observe is a stiffening, rather than 
a softening,^of the resonance. We conclude therefore that the inertial term from linear 


momentum 


Ref. j^. The former term is a direct product of the breathing Fermi surface model causing 
conductivity-like Gilbert damping jlR. 17|. In the conductivity-like model the value of ob¬ 
served r, 0.1-0.5 ps, represent the average of Bloch state lifetime tb near the Fermi surface 
for Py and Co. 


, y, la] dominates over the inertial term from angular momentum proposed 


m 




FIG. 3. Ferromagnetic resonance half-power linewidths po^Ri /2 a function of a;/27r for (a) Py 
and (b) Co, between 4 to 345 GHz. Dashed lines are linear extrapolations of the low-frequency 
FMR linewidths; solid curves are fits to Eq. (3). Each point with error bar represents a summary 
of six measurements. 
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Note that Gilbert damping coefficient a, calculated by the breathing Fermi surface mode 


should also be frequency-dependent. A relation oi a = a^/{1 + rl) can be obtained 


Q 


16 


with Go the zero-frequency damping. Fig. 3 shows the linewidths of both low-frequency 
and high-frequency ranges for Py and Co. Here points with errorbars (one a) are used for 
frequencies > 100 GHz to reflect the scattering of six repeats. The dashed lines are low- 
frequency linear extrapolations from oq. Compared with extrapolated linewidths from the 
low-frequency AiPi/ 2 (cu), we observe reduced linewidths at high frequencies (115-345 GHz), 
icit prediction has been made for the effect of rotational inertia on the linewidth in 


No exT 
Ref. t 


12 |. but 


terms in Ref. 


he observed behavior matches well with the prediction of linear inertial 


16l-ll8|. The solid curves are fits to the form: 

2q;oi^ 1 


PoAhfi/2 — /ioAi/o 


7 1 - 1 - 


( 3 ) 


taking the ccq from the low-frequency linewidths. The fitted r are listed in Table 1. The 
values of r extracted from /ioAFri /2 are close but slightly smaller than from /loFTres- We 
note here that the linewidth measurements are less precise than /ighfres and more sensitive 
than resonance fields to various sources of inhomogeneities in the samples. 

At room temperature, the relaxation rate of Bloch states (1/tb) is determined by the 
electron scattering with phonons and impurities. In this sens e, tr is similar to the remagneti¬ 
zation time Tr in the ultrafast demagnetization experiments 191421|. where optically excited 
electrons also relax through electron-phonon and electron-impurity interactions. In the limit 
of zero laser fluence, nonlinear effects due to high occupation number of excited states are 
reduced}^. 22|. The zero-fluence tr has been reported to be 0.2-0.25 ps for Py 10-30 nm 22] 
and ~ 0.4 ps for Go 15 nm 2l(|, close to the value of r in this work. 

At high frequencies the (nonmagnetic) skin depths 5s in the ferromagnetic films become 


smaller and the enhanced eddy current effect may influence the resonance field 23j . The res¬ 
onance field will be enhanced by ~ fioMs/S^k^, where feg is the lowest-energy wavenumber 
determined by the surface anisotropy. Because 5^ oc l/cu, the resonance field enhancement 
is proportional to and may influence the quadratic term in Eq. (2). Our calculations 
show that this term is negligible, about 0.4 mT for Py 30 nm and 0.09 mT for Go 30 nm, 
compared with the observed effects up to 80 mT (See the Supplemental Information for 
details). 

We do not believe that interfacial effects, including spin pumping, play an important role 





















in the observed high-frequency behavior. Both Gilbert damping ao and the two inertial 
dynamics lifetimes r^fioHres) and r(/ioAhfi/ 2 ) in Table 1 show little thickness dependence 
for either Py or Co, indicating that bulk relaxation is dominant. The weak thickness de¬ 
pendence of (To is consistent with the very low spin pumping effect of Py/Ta identihed in 


Ref. 


24l |. in any case negligible for 30 nm dims and without quadratic frequency depen¬ 


dence. Theoretical predictions of resonance shifts from imaginary spin mixing conductance 


are three orders of magnitude lower than observed here 


25l |. Only r(/ioh^res) measured from 


resonance shifts for 6 nm dims, not matched in T{fiQAHi/ 2 ) measured from linewidths, dider 
signidcantly. This enhancement may be structural in origin. 

A technological implication of our results is that edective deld requirements for proces¬ 
sional switching will be reduced as switching times in magnetic storage decrease into the 
few-picosecond range. In this sense, the inertial dynamics ease ultrafast switching, if the 
behaviors of Py and Co are representative of other metallic ferromagnets. The edective deld 
reduction, here up to 80 mT in Co, is not small in an absolute sense, and might according 
to Eq. (2) be enhanced signidcantly in ferromagnets with higher Gilbert damping. On the 
other hand, prior switching experiments with high-held relativistic electron bunches seem 
to indicate that nonlinear damping increases edective deld requirements by a rather larger 
amount for large-angle dynamics in CoCrPt 8j, underscoring the utility of HF-FMR to iden¬ 
tify the inertial edect. 

In summary, we identify a novel term to magnetization dynamics in the ferromagnetic 
metal dims Py and Co which is quadratic in frequency and becomes significant above 100 
GHz. The term stidens the frequency, aiding applied delds in driving ultrafast magnetization 
motion. The behavior is best explained W ^namics retarded through a finite Bloch-state 
relaxation time tb as proposed in Refs. 3|, y, l6|. Extracted relaxation times are 0.1-0.2 ps 
for Py and 0.2-0.4 ps for Co, close to the remagnetization times measured in optical pump- 
probe demagnetization experiments. Our dndings extend understanding of magnetization 
dynamics at picosecond time scales and may open up new possibilities for high-speed inertial 
switching in ferromagnetic materials, previously demonstrated only in antiferromagnets 10|. 
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